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ABSTRACT 

We have measured a time delay of 26^ 4 days and a magnification ratio of 1.52 ±0.05 
in the strong radio gravitational lens PKS 1830—211. The observations were made 
over the 18 month period from 1997 January to 1998 July with the Australia Telescope 
Compact Array at 8.6 GHz, and have shown that the source started a large flux density 
outburst around 1997 June. 
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1. Introduction 



A precise measurement of the time delay between the images in a gravitational lens is essential 
if it is to fulfill its potential as a tool for estimating the Hubble constant. The flat-spectrum radio 
source PKS 1830—211 ( Rao fc Subrahmanyan 1988j ) was found to be an Einstein ring/gravitational 
lens, comprised of two compact, flat-spectrum components located on opposite sides of a 1 arcsec 
diameter ring ( Jauncey et al. 199 1| ). It is by far the strongest radio lens yet found, being ~ 10 Jy 
at 4.8 GHz. Observations of molecular absorption at centimeter and millimeter wavelengths have 
revealed two intervening galaxies at z = 0.19 (Lovell et al. 1996a) and z = 0.89 ( [Wiklind fc 



Combes 1996| ). VLBI observations of the two compact lensed components have revealed striking 
structural differences on milliarcsecond scales ( Garrett et al. 1996| ) suggesting that both galaxies 
may be involved in the lensing (Lovell et al. 1996a). PKS 1830—211 varies dramatically at radio 
wavelengths flLovell et al. 1996b ) and VLBI observations show that the variations are confined to 



the two compact components (King 1994), making it an excellent candidate for relative time delay 
measurements. 

This paper describes our measurements with the 6 x 22 m Australia Telescope Compact 
Array (ATCA) to determine the lensing time delay in PKS 1830—211 through a correlation of the 
flux density light curves of its compact components. Section 2 describes the development of our 
observing strategy and data reduction, section 3 describes the analysis of our data to obtain a 
time delay and magnification ratio. In section 4 we describe our error analysis and in section 5 we 
discuss the implications of our measurements on modeling the PKS 1830—211 lensing system. 



2. Observations and Data Reduction 

The presence of a single peak in the 8.4 GHz total flux density light curve in early 1992 
( Lovell et al. 1996b| ) indicates that the lensing time delay is less than a few months (a time delay 
that was much longer would result in a broader or double-peaked light curve). Consequently, in 
order to sufficiently time resolve these variations in the compact components it was necessary to 
measure the flux densities of the two components separately on timescales shorter than both the 
total flux density variations and the lensing time delay, previously estimated to be 44 ±9 days (van 
Ommen et al. 1995). We commenced our ATCA monitoring observations in 1995 August, making 
observations at approximately monthly intervals. The two compact and variable components 
can be measured separately, as explained below, on the 6 km East- West array of the ATCA at 
8.6 GHz if observations are restricted to between 15 and 19 hours LST when the ATCA provides 
sub-arcsecond spatial resolution along the position angle of the compact components. The strong, 
compact, flat-spectrum source PKS 1921-293 was observed in each session to provide the relative 
gain calibration of the individual antennas, while overall flux density calibration was determined 
through observations of the ATCA primary flux density calibrator PKS 1934-638. The total flux 
density variations of PKS 1830—211 will be presented in a future paper (Lovell et al. 1998, in 
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preparation). During the period 1996 July 26 through 1996 August 8, a more closely spaced series 
of observations was undertaken in order to gain an understanding of the shorter-term variability. 
Analysis of this 1996 July- August data showed more rapid variability than previously thought, so 
from 1997 January, one observation every three to six days was made to better sample the light 
curves. In this paper we restrict our analysis to the data from 1996 July onwards, although the 
data from the period 1996 December 20 through 1997 January 10 were also excluded as the source 
was too close to the sun. 

As the two compact components of PKS 1830—211 are not fully resolved from the ring at the 
~ 0.9 arcsecond resolution of the ATCA at 8.6 GHz, we have adopted a model fitting approach to 
determine the flux densities of the individual components. We chose a simple model consisting of 
two circular Gaussian components, each of which also contains a contribution from the non- varying 
ring. The component separation was set initially to 0.95 arcseconds at a position angle of 46 



degrees, as determined from VLBI observations (King 1994). These values were allowed to change 
slightly during model fitting, but were found to deviate by no more than 3 degrees in position 
angle and 12 mas in separation. All data were edited and calibrated using AIPS before being 
exported to Difmap ( |Shepherd 1997] ) for model fitting. 

To verify that this model successfully describes the source, we analyzed an edited section 
from a full ATCA synthesis observation taken on 1994 December 19 when PKS 1830—211 had a 
total flux density of 7.2 Jy. We found that the individual fitted Gaussian component flux densities 
differed by no more than 70 mJy from those measured from the full synthesis data. As a further 
check, simulated ATCA data sets were generated and analyzed using the above procedures. For 
all of these the modeled component flux density changes were accurate to within 1%, and the total 
flux density in the model accounted for 99.8% of the total simulated flux density. This clearly 
demonstrates the robustness and reliability of our model fitting procedures in determining small 
changes in each of the compact components of PKS 1830—211. 

Component flux density errors were estimated at each epoch using the method described by 
Tzioumis et al. (1989) and were found to be typically ~ 50 mJy. The NE and SW component flux 
density light curves obtained from the fitted 8.6 GHz ATCA data since 1996 January are shown 
in Figure 1(a). We have detected clear changes in the flux densities of the two components, with 
a roughly linear increase seen to begin sometime between MJD 50530 and MJD 50670. There is 
a clear difference in gradient between the two light curves, which is a result of different image 
magnifications. The figure also shows some significant differences in the light curves, the most 
prominent being a distinct 'bump' of 400 mJy in the NE component light curve near MJD 50700 
that is seen to occur significantly later in the SW light curve, indicating a time delay of 20 to 30 
days. In the following section, we discuss a detailed analysis of the data to obtain a quantitative 
estimate of the time delay and magnification ratio and their uncertainties. 
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3. Analysis 

To estimate the time delay it is important to consider how the variability of the lensed source 
is detected by the observer. In the following analysis, the subscripts 1 and 2 refer to the NE and 
SW components of PKS 1830—211 respectively. We describe the time variation of the NE and 
SW compact components of the lensed image as Si(t) and S 2 (t) respectively. The SW component 
light curve is expected to be identical to the NE light curve but magnified by a different amount 
and shifted in time, thus: 

S 2 (t) = -Sx{t + Ar) (1) 

where At is the time delay and \i is the relative magnification ratio. Also lensed is an extended 
component that forms the Einstein ring which contributes a constant flux density S^ng- ^he total 
observed flux density can therefore be described as 

S m (t) = S^t) + S 2 (t) + S ring . (2) 

The Einstein ring contributes a large proportion, up to 50% of the total flux density at cm 
wavelengths ( pauncey et al. 1991 ). Therefore the effect of the ring flux density must be considered 



before attempting to determine a lensing time delay. In our ATCA observations, neither Si(f) nor 
S 2 (t) are measured directly, instead the measured quantities are 

S ml (t) = S^t) + S cl (3) 

and 

S m2 {t) = S 2 (t) + S c2 = -5i(t + At) + S c2 (4) 

A* 

where the quantities S c \ y2 are constant but unknown and 5" r i n g = S& + S c2 . It is unlikely that S c \ 
and S c2 are equal. Therefore it is also necessary to solve for d, the difference between S c i and S c2 . 
The two light curves may now be compared using trial values of /i, At and d so that a solution 
for the relative magnification ratio and the time delay can be obtained. 



3.1. Dispersion Analysis 

A correlation of the two light curves was undertaken using the dispersion analysis method 
introduced by Pelt et al. (1994, 1996) to analyze the component light curves of the gravitational 
lens 0957+561. The dispersion method has been shown to be successful for 0957+561, which has 
a relatively large time delay of 417 ± 3 days ( Kundic et al. 1997 ). We chose this method because 
it avoids any interpolation between data points. Methods that are based on interpolation can lead 
to erroneous results as they put equal weight on assumed and measured data. 



Following Pelt et al, for every test value of fj,, At and d, two light-curve datasets ctj and 
bi (i = 1, . . . , N where N is the number of observations) are obtained. Together they form a 
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combined light-curve where ctj contains the NE variable component data and bi contains the SW 
variable component data modified in flux density and time by the test values for fi, Ar and d. 
The dispersion, D 2 , of this combined light curve is calculated from the weighted sum of squared 
differences between nearby ai,bi pairs over the entire curve: 



where 



^ {AT ^ = ^3hzhl (5) 

1 if \ti — tj\ < 5, 

K = { , 2N -i («) 

(l+(^7^) ) if\ti-tj\>5. 



W ZJ = (7) 



and Wi = 1/ef where is the standard error in data point i. The weighting factor V/j 
ensures that data pairs are only given significant weighting if they are less than ~ 5 days apart 
and is modified from the Vij of Pelt et al. to avoid strong weighting on data points with near-zero 
separation. These points add a bias towards a zero delay for relatively short time delays, as we 
expect to find in PKS 1830—211. The values of Ar, n and d for which D 2 is a minimum are thus 
the best estimate of the time delay, magnification ratio and constant flux density offset difference. 

Before a correlation of the light curves can be made, an estimate of 5, effectively the longest 
timescale over which the source is believed not to show significant variability, must be determined. 
An inspection of the light curves in Figure 1, shows significant changes in flux density on timescales 
of the order of 10 days. We note that V' drops quite slowly, reaching 0.5 at a point separation 
of 1.55. Therefore, we chose a value of 5 = 5, slightly less than the estimated timescale, so that 
points with separations greater than ~ 8 days did not unduly influence the solution. 

A single, unambiguous solution was found for 5 = 5 days at Ar = —23 days (i.e. NE 
component leading), = 1.52 and d = —0.62 Jy. We also investigated how the solution behaved 
for a range of 5 and found that for 5 < 7 days, Ar stayed within ±4 days, n stayed within ±0.02 
and d stayed within ±0.06 Jy of the 5 = 5 days solution. For 5 > 7 days, up to 6 = 14 days (the 
largest value we trialed), the solutions for Ar became shorter, ranging from —12 to —17 days, but 
solutions for (i and d changed little. A visual inspection of the combined light curves with these 
solutions applied clearly shows a poorer correlation than seen for 5 < 7. 

It is likely that this sudden change in the delay solution is due to an effective smoothing of 
the features that help register the light curves in delay. However, \i and d are affected to a lesser 
degree as the overall smooth, rising light curves constrain these parameters well, and the short 
timescale changes have little influence. 
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4. Error Analysis 

We carried out Monte Carlo simulations to estimate the errors in our derivation of the time 
delay and magnification ratio. Five hundred light curve pairs were created with data points 
at the same epochs of the actual observations. The flux density of each point was calculated 
using a random number generator but constrained such that each simulated point was within the 
uncertainty of the original measurement. An additional constraint forced the sum of a pair of 
same-epoch measurements to be equal to the total measured flux density within the uncertainty 
of that observation. 

Each of the 500 light curve pairs were then correlated to obtain estimates for At, and 
d when 5 = 5 days. We found that the correlated values of \i and d were distributed in a 
Gaussian-like manner about the values obtained from our observations. We have used these 
distributions to estimate errors for n and d by measuring the width of the distributions at half 
their peak. We thus estimate that the magnification ratio, /i, is 1.52 db 0.05 and the flux density 
offset, d, is —0.62 ±0.05 Jy. The distribution of At from the simulations however is not distributed 
in a Gaussian-like fashion. Instead, a broad distribution of solutions between -12 and -30 days 
was obtained with peaks near -15, -23 and -28 days. 

Eighteen years of monitoring the gravitational lens 0957+561 have demonstrated that 
non-varying or smooth light curves contributes little information to a time delay solution; in fact 
such data can lead to ambiguous solutions. The detection of strong features in both light curves is 
required to obtain a reliable time delay measurement flKundic et al. 1997| ). Bearing this in mind, 
we have also analyzed the subset of data surrounding the 'bump' seen near MJD 50700 where the 
most prominent and best sampled feature occurs (Figure 2(a)), specifically to avoid ambiguities in 
the time delay estimate. 

We found a solution at At = — 26 days, \x = 1.68 and d = — 1.17 Jy for S = 5 (Figure 2(b)). 
As expected, Monte Carlo simulations show the magnification ratio and flux density offset to be 
less well constrained here than with the full dataset. We have therefore adopted the values for \x 
and d as determined from the full dataset as our best estimates. However, the simulations show 
the time delay for this subset to be well constrained without ambiguities with a Gaussian-like 
distribution of solutions around -26 days. By measuring the width of this distribution at half its 
peak, we have estimated an uncertainty and obtain a time delay of —26^4 days. 

We are continuing our ATCA monitoring observations. When the current trend of increasing 
flux density at 8.6 GHz ends, the change in gradient may provide even better constraints on the 
time delay. 
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5. Discussion 

Out ATCA monitoring observations have allowed us to measure the time delay and 
magnification ratio of the gravitational lens PKS 1830—211, providing new constraints for lensing 
models of this system. 

Our time delay measurement is not in agreement with the —44 ± 9 day measurement of van 
Ommen et al. (1995) from VLA monitoring observations at 8.6 and 15 GHz. We believe van 
Ommen et al. did not correctly account for the contribution of the Einstein ring flux density when 
calculating the magnification ratio. The difference between the total flux density of the source 
and the sum of the compact component flux densities fitted to the data, i.e. the contribution from 
the ring, should be constant. However in van Ommen et aVs analysis the inferred Einstein ring 
contribution depends on the observing frequency and array configuration. We feel confident that 
our analysis of the ATCA data has correctly accounted for the constant Einstein ring flux density. 

Millimeter observations, in early 1996, by Prye, Welch, & Broadhurst (1997), which clearly 
separate the core components of the lens and show no extended jet structure on ~ 0.5 arcsec 
scales, provide an estimate for fx of 1.14 ± 0.06. Also, Wiklind & Combes (1998) measured the 
magnification ratio in the millimeter band on three occasions over 11 months and found that it 
changed between 1.0 and 1.8 with an uncertainty in each case of 0.1. We note that PKS 1830—211 
has shown strong variability at millimeter wavelengths in the past ( [Tornikoski et al. 1996| ) 
and that Wiklind & Combes observed a decrease in total flux density from 1.5 to 0.93 Jy over 
their observations. Such variability could influence the measurement of the magnification ratio 
at a single epoch. Also, the likelihood of sub-components in the compact structure at these 
wavelengths, similar to those seen at 43 GHz by Garrett et al. (1998), should be investigated 
as this too has an impact on the magnification ratio. More detailed monitoring of the compact 
components at millimeter wavelengths can provide a test of our magnification ratio measurement 
if the time delay effect is accounted for. 

It is not yet possible to obtain an estimate for Ho from the PKS 1830—211 lensing system as 
the redshift of the lensed AGN is not known, and the lensing galaxies are not well parameterized. 
However, if we assume that the intervening galaxy at z = 0.19 has a negligible gravitational affect 
on the system, then we can apply our time delay measurement to the model of Nair et al. (1993) 
for a lensing galaxy at z = 0.89 to put constraints on the redshift of the background source, z s . If 
Ho is assumed to be between 50 and 100 kms _1 Mpc _1 and qo = 1/2, then a time delay of —26^ 
days implies z s > 1.4, with a value of z s ~ 2 for Ho = 75 kms _1 Mpc _1 . 

The Australia Telescope is funded by the Commonwealth Government for operation as a 
national facility managed by CSIRO. J.E.J.L. acknowledges the receipt of a JSPS Fellowship. 
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Fig. 1. — (a) The 8.6 GHz light curve data for both components of PKS 1830—211. The measured 
NE component flux densities are represented by the open triangular symbols and the measured 
SW component flux densities are represented by the closed circular symbols, (b) The NE and SW 
component light curves after applying our solutions from the dispersion analysis (At = —23 days, 
H= 1.52 and d = -0.62 Jy). 
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Fig. 2. — (a) The 8.6 GHz light curve data for both components of PKS 1830—211 for the time 
period surrounding the 'bump' described in Section ^. The measured NE component flux densities 
are represented by the open triangular symbols and the measured SW component flux densities 
are represented by the closed circular symbols, (b) The NE and SW component light curves after 
applying our solutions from the dispersion analysis (At = —26 days, \x = 1.68 and d = —1.17 Jy). 



